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Abstract 

The Grobner walk is an algorithm for conversion between Grobner bases for different 
term orders. It is based on the polyhedral geometry of the Grobner fan and involves 
tracking a line between cones representing the initial and target term order. An im- 
portant parameter is explicit numerical perturbation of this line. This usually involves 
both time and space demanding arithmetic of integers much larger than the input num- 
bers. In this paper we show how the explicit line may be replaced by a formal line 
using Robbiano's characterization of group orders on Q". This gives rise to the generic 
Grobner walk involving only Grobner basis conversion over facets and computations with 
marked polynomials. The infinite precision integer arithmetic is replaced by term or- 
der comparisons between (small) integral vectors. This makes it possible to compute 
with infinitesimal numbers and perturbations in a consistent way without introducing 
unnecessary long integers. The proposed technique is closely related to the lexicographic 
(symbolic) perturbation method used in optimization and computational geometry. We 
report on computations with toric ideals, where a version of our algorithm in certain 
cases computes test sets for hard integer knapsack problems significantly faster than the 
Buchberger algorithm. 

1 Introduction 

Let R = k[xi, . . . , Xn] denote the polynomial ring in n variables over a field k. Grobner basis 
computations in R tend to be very expensive for certain term orders (like the lexicographic 
order) . Therefore it often pays to compute Grobner bases for "easier" term orders and convert 
them into Grobner bases for the desired term order. For zero- dimensional ideals this can be 
accomplished by the FGLM-algorithm [Hj. For general ideals the Grobner walk algorithm [H] 
can be applied. 

Let -<i and -<2 be term orders on R. The usual Grobner walk proceeds from the reduced 
Grobner basis G for / over -<i by tracking a line u;(t) = (1 — t)ujQ + tTQ, < t < 1, where 
ti^o and To are vectors in the respective Grobner cones C^-^{I) and C^^i^) of I. At t = the 
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Grobner basis is known. The line uj(t) is tracked through the Grobner fan of I and Grobner 
bases are computed at common faces of successive Grobner cones. At t = 1 we reach the 
reduced Grobner basis for / over -<2. 

The efficiency of the Grobner walk rests on clever choices of c^o and tq. A choice of uq and tq 
on low dimensional faces of Grobner cones may lead to very heavy Grobner basis calculations 
along uj{t). Often (but not always) it pays to choose Uq and Tq generically inside C_<j(/) and 
C^^{I) and ensure that u!{t) only intersects common faces of low codimension on its way to 
the target term order -<2. 

The initial reduced Grobner basis G over -<i makes it possible to compute an interior 
point in C_<j(J). Computing an interior point in the target cone C^^{I) is considerably more 
difficult, since we do not know the reduced Grobner basis over -<2 in advance. Tran jTH] 
approached this problem using general degree bounds on polynomials in Grobner bases. The 
general degree bounds in Tran's approach may lead to integral weight vectors with 10, 000- 
digit entries in representing a lexicographic interior point in the case of polynomials of degree 
10 in 10 variables. 

In this paper we give an algorithm where the line uj{t) is replaced by a (formal) line VL{t) 
between suitably chosen perturbations given by and -<2 and /. It turns out that the 
numerical dependence on / disappears in our algorithm and that VL[t) may be viewed as a 
line which can be used for all ideals in the Grobner walk from to -<2- The formal line has 
the property that its initial and target points are always in the interior of the Grobner cones. 
Furthermore the common faces that Vt{t) intersect are all facets. 

In the classical Buchberger algorithm j2] for computing Grobner bases one only computes 
with term orders and initial terms of polynomials. Tracking Vt{t) gives a "Buchberger- like" 
Grobner walk algorithm, where one only needs to compute with marked polynomials and term 
orders. On smaller examples the algorithm can easily be carried out by hand (cf. 

We have observed some interesting experimental results using a version of the generic walk 
tailored to lattice ideals [12]. When the generic walk is applied in computing full test sets 
for feasibility of the hard integer knapsacks from the natural initial and target vectors 
are rather close in the Grobner fan. This leads to very fast computations of test sets. These 
examples with polynomials of high degree in many variables seem out of reach for the classical 
Grobner walk. We report on computational experiments in the last section of this paper. 

An understanding of our algorithm requires a firm grip on the usual Grobner walk algo- 
rithm. Therefore ^and ^recalls and proves fundamental results for the usual Grobner walk 
using which we transition to the generic Grobner walk in 21 

The basic technique we propose to avoid explicit perturbation is not quite new. The key 
idea of implicit (symbolic) perturbation was proposed by Charnes in 1952 to make Dantzig's 
simplex method for linear programming finite. The method is now known as the lexicographic 
perturbation method, see |31 page 34], and used by many reliable implementations of the 
simplex method. In computational geometry, similar symbolic perturbation schemes are used 
to treat input data points in M" as if they were in general position, see [HI page 14]. 
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2 Preliminaries 

In this section we recall the basics of convex polyhedral cones. We emphasize a crucial result 
from the theory of group orders (Lemma 12.11) and recall the construction of the (restricted) 
Grobner fan. 

2.1 Cones and fans 

A convex polyhedral cone is a set 

C{Vi, ...,Vr)= R>oVi + ■■■+ R>oVr C W 

where Vi, . . . ,Vr G M*^. In the following a cone will refer to a convex polyhedral cone. The 
dual of a cone C C is 

= {cj G I {uj, v) > 0, for every v e C}. 

The dual of a cone is a cone and the intersection of two cones is a cone. The dimension 
of a cone is the dimension of the linear subspace it spans. For a vector u G M" we let 
M"^ = {x G I (m, x) = 0}. A face F C C of a cone C is a subset F = w^nC, where u G C^. 
Faces of codimension one in C are called facets. 

A collection JF of cones and their faces is called a fan if for every Ci,C2 G we have 
Ci n 6*2 G JF and Ci fl C2 is a common face of Ci and C2. 

2.2 Rational group orders on Q"^ 

Let [A, +) be an abelian group. Recall that a group order -< on A is a total order -< on A 
such that 

X ^ y =^ X + z ^ y + z 

for every x,y,z & A. 

Let u = {ui, . . . , Un) C Q"' be a Q- vector space basis for Q". Then we get a group order 
on Q"' given by u v if and only if 

{{Ui, U),..., {Un, U)) <iex {{uJi, v) , . . . , (u;„, v)) , 

where <iex refers to the lexicographic order on Q". We call such a group order rational. To 
describe arbitrary group orders on similarly, one needs a more general setup including real 
vectors (see fE])- To ease the exposition we will restrict ourselves to rational group orders. 
A group order refers to a rational group order in the following. For a rational e > we put 

Ue = UJi + eLU2 -\ h e"'~^UJn. 



The following well known lemma plays a key role in the generic Grobner walk. 
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Lemma 2.1 Let uo = (cui, . . . C Q" be a Q-basis. Suppose that F C Q" is a finite set of 
non-zero vectors with v for v & F . Then there exists < 5 G Q such that {u^, v) > for 
every v E F and e G Q with < e < 5. 

Proof. We prove this by induction on n. The case n = 1 is clear. For n > 1 we may find 
< (5o G Q such that 

{Un^l + eUn, f ) > 

for every v E F with {un-i,v) > and e G Q with < e < 6o. Therefore -<uj' v for 
Lj' = {lui, . . . , Un-2, + ^<^n) for cvcry f G F if < e < (5o. 

By induction there exists < 5i G Q such that (ci;^, f ) > for every v E F and e G Q 
with < e < 5i. Putting 5 = min((5o, Si) we get (a;^, f ) > for every v E F and e G Q with 
0<e<6. □ 

A group order -< on Q" is called a term order if -< f for every G N". This is equivalent 
to ^ Cj where denotes the i-th canonical basis vector for i = 1, . . . ,n. As a consequence 
of Lemma 12.11 we get the following corollary. 

Corollary 2.2 Let F C Q" be a finite set of positive vectors for the group order -< i.e. v >~ Q 
for every v E F . Then there exists uj E Q" such that 

{u, v) > 

for every v E F. If ^ is a term order, we may assume that uj has positive coordinates. 
2.3 The Grobner fan 



Let R = k[xi, . . . , Xn] denote the ring of polynomials in n variables over a field k. It is 
convenient to view R as the semigroup ring A;[N'^]. We briefiy recall the construction of the 
(restricted) Grobner fan (cf. ^2]) for an arbitrary ideal in R. 

Fix a group order -< on Q*^. For a polynomial / = ^^gj^n a^x'" E Rwe let supp(/) = {v E 
N"' \ tty 0} and in_.(/) = a„a;", where u = max^ supp(/). For a vector a; G M" we let in^(/) 
denote the sum of terms a^x^ in / maximizing the u-weight {u, v). We call / uj-homogeneous 
if / = in^(/). To an ideal I C R we associate the ideals in^(J) = (in^(/) \ f E I \ {0}) and 
in^(J) = (in^(/) | f E I). These ideals may be viewed as deformations of the original ideal J. 
The initial ideal in_.(J) is generated by monomials. This does not hold for in^(J) in general 
(unless u is chosen generically) . 

Now define 

dM) = {u- I u' E supp(/) \ {u}} C Z^ 
where QuX^ = in^(/). For a finite set F C R of polynomials we let 

and 

C^(F)=C(9^(F)rnM^o 

= {ujE M^o I ^) > 0, G 9^(/), / G F}. 
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Notice that dimC^(F) = n by Corollary 12.21 and that 

C^{F) = {ue Rl, I in^(in^(/)) = in^(/) for every / G F}. 

A Grobner basis for I over ^ is a finite set of polynomials G = {gi, . . . , Qr} C J such that 

(in^(5fi), . . . , in^igr)) = in-<(/). 

The Grobner basis G is called minimal if none of gi, . . . ,gr can be left out and reduced if 
the coefficient of m^{gi) is 1 and in^(5fj) does not divide any of the terms in gj for i ^ j 
and i, j = 1, . . . ,r. A reduced Grobner basis is uniquely determined. Minimal Grobner bases 
exist for arbitrary group orders. However, Grobner bases over arbitrary group orders do not 
necessarily generate the ideal (as opposed to Grobner bases over term orders). Similarly, the 
reduced Grobner basis is only guaranteed to exist for term orders. 

To define the Grobner fan we now specialize to the case where -< is a term order. The 
Grobner cone G^{I) of an ideal / over -< is defined as C^(G), where G is the reduced Grobner 
basis of I over -<. The Grobner fan of / is defined as the set of cones G^{I) along with their 
faces, where -< runs through all term orders. This is a finite collection of cones jl7| Theorem 
1.2] and one may prove that it is a fan (Propositions 2.3 and 2.4 in ^\ give a proof assuming 
non- negative weight vectors). The following proposition shows that G^{I) is the largest cone 
among G^{G), where G is a Grobner basis for I over -<. 

Proposition 2.3 Let G be a (not necessarily reduced) Grobner basis for I over -<. Then 

G4G)CG4I). 

Proof. Transforming G into a minimal Grobner basis G' by omitting certain polynomials in 
G clearly leads to an inclusion G^{G) C G^{G'). Transforming G' into the reduced Grobner 
basis proceeds by a sequence of reduction steps: suppose that /«, fj G G' and that a term 
in fj is divisible by in^(/j). Then fj is replaced by /j = fj — {x'" / in^(/j))/j. This reduction 
may introduce "new" monomials which are not present in fj. More precisely if w G supp(/j), 
then w G supp(/j) oi w = v — u + u' , where ttuX'^ = in^(/j) and u' G supp(/j). In the 
latter case we get w' — w = {w' — v) + {u — u'), where a^/x"' = m^{fj). Let G" denote the 
Grobner basis obtained by replacing fj with /j. Then G{d^{G')) 3 G{d^{G")) and thereby 
C^{G') C G^{G"). Since the reduced Grobner basis is obtained using a finite number of these 
reduction steps we have proved the inclusion. □ 

For a specific term order one may have infinitely many cones given by different minimal 
Grobner bases. As an example consider the ideal I = {x, y) C k[x, y]. If n is a positive natural 
number then Gn = {x — y", y} is a minimal Grobner basis for / over the lexicographic order 
-< with X >~ y. In this case 

GM)^C^iG^)^G^iG2)^■■■ . 

3 The Grobner walk 

We outline the basic idea of the Grobner walk [H] and give a new lifting step using reduction 
modulo the known Grobner basis. 
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Let -<i and ^2 be term orders and / an ideal in R. Suppose that we know the reduced 
Grobner basis G for I over ^1. If 

u;eC^,(/)nC^,(J) 

hes on the common face of the two Grobner cones, then G^^ = {ini^{g) \ g G G} is the reduced 
Grobner basis for in(^(/) over -<. Now a "hfting" of G^j to a Grobner basis for / over -<2 is 
required. The procedure for this is based on Proposition IH.2I below. It involves a Grobner 
basis computation for m^{I) over -<2. The point is that if F = C^^(/) fl C^^{I) is a high 
dimensional face (like a facet) and u is in the relative interior of F, the ideal in^(/) is close 
to a monomial ideal and this Grobner basis computation becomes very easy. 

Given a term order -< and a vector uj G M>o define the new term order by u v if 
and only if (m, uj) < (f , uo) or (m, uj) = (f , uj) and u -< v. We record the following well known 
lemma. 

Lemma 3.1 Proposition 1.8] Let I C R be any ideal and uj G M>g. Then 

in^(in<^(/)) = in^^(/). 

The lifting step (Proposition IH.2r ii) below) in the following proposition is different from 
the lifting step in the usual Grobner walk jTTl Subroutine 3.7]. 

Proposition 3.2 Let I C R be an ideal and -<i, -<2 term orders on R. Suppose that G is the 
reduced Grobner basis for I over ^1. If uj ^ G^^{I) fl G^,^{I), then 

(i) The reduced Grobner basis for in^(/) over ^1 is G^i = {in^(5f) \ g ^ G}. 

(a) If H is the reduced Grobner basis for in^^^I) over -<2, then 

{/ -fife H} 

is a minimal Grobner basis for I over <2uj- Here f'^ is the unique remainder obtained 
by dividing f modulo G. 

(Hi) The reduced Grobner basis for I over -<2i^ coincides with the reduced Grobner basis for 
I over -<2. 

Proof. Given a term order -< and a vector uj G G^{I), the reduced Grobner bases for / over 
-< and agree. This proves (iii) and (i) taking Lemma 13.11 into consideration. Suppose 
that / is an u;-homogeneous polynomial (cf. §2.3|1 in int^(J). Using the division algorithm in 
computing the unique remainder f^, we keep reducing terms with the same cu-weight as the 
terms in / = in^(/). Since m^-^{m^{g)) = in^-^(5f) for G G and /"^^ = 0, we see that all 
terms in will have cj-weight strictly less than the terms in /. Therefore 

in^(/) = in^(/-/^). 
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Now suppose that {fi, ■ ■ ■ , fs} is the reduced Grobner basis for m^{I) over -<2. In particular 
we get that fi is cu-homogeneous for i = 1, . . . , s. Then 

in^2<^(/) = in^2(in^(/)) = (in^2(/i), . . . , in^2(/,)) 
= (in^2(in<^(/i)), . . . , in^2(in<^(/,))) 
= (in^,(in^(/i - f^)), . . . ,in^,(in^(/, - /f ))) 
= (in^,JA-/n,...,m^.J/s-/f)). 

This proves that {/i — /f , . . . , fs — f^} C J is a (minimal) Grobner basis for I over -<2uj- □ 

Proposition 13.21 may be turned into a Grobner basis conversion algorithm as shown in the 
following section. 

3.1 Conversion along a line 

A natural approach to Grobner basis conversion is to trace the line between vectors in different 
Grobner cones and update Grobner bases successively using Proposition 13.21 This process is 
called the Grobner walk [Sj. A good reference for this procedure is [HI §4], which inspired the 
following. We sketch the first step of the Grobner walk. The succeeding steps of the Grobner 
walk are similar. Suppose that ujq G C^-^{I),Tq G C^^{I) and that G is the reduced Grobner 
basis for I over ^i. Here and -<2 are rational term orders (cf. ^2.2j) given by Q-bases 
uj = {uji, . . . , ujn) and r = (ti, . . . , r„) respectively. Then we consider the line 

u;{t) = (1 - t)iUo + tTo, 0<t<l 

in the Grobner fan of / from uq to tq. Initially we know the reduced Grobner basis at a;(0) = uq 
(being G). Consider the "last" u' = uj{t') in C^-^(J) = C^^{G). To be more precise t' satisfies 

1. < t' < 1 

2. uj{t) G C^i(/) for t G [0,t'] and uj{t' + e) ^ C^^{I) for every e > 0. 

If no such t' exists then G is the reduced Grobner basis over -<2. If t' exists oj{t') is on a 
proper face of C^^{I) and v G d{G) exists with {uj{t' + e),t>) < for e > 0. This implies that 
{to,v) < {ujq,v) and hence (ro,f) < 0. 

This indicates the procedure for finding t' given G. For v G d{G) satisfying {to,v) < we 
solve {uj{t), v) = for t giving 

{uJo,v) - {To,vy 

Then t' is the minimal among these ty. In this case u' = uj{t') lies on a proper face F of C^^{I) 
and clearly 

Now we use -<2u' the term order -<2 in Proposition 13.21 The point is that we only need 
the target term order -<2 to compute a Grobner basis for in^/(J) (not the notational beast 
-<2uj')- The reason for this is that the Buchberger algorithm in this case solely works with 
cj'-homogeneous polynomials and ties are broken with -<2. 
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To prove that we actually enter a new Grobner cone we need to show that t' > (cf. [H], 
§5, (5.3) Lemma). In the initial step it may happen that t' = 0. But if this is the case we 
may assume (in the following step of the Grobner walk) that G is the reduced Grobner basis 
over -<= -<2ui'- Since tq G C^^i^) is non-zero we may use tq as the first vector in a Q-basis 
representing -<2. In this case assume that t' = 0. This means that we can find v G d{G) with 
(a;', v) = and (tq, f ) < contradicting that G is a Grobner basis over -<. 

We have outlined the procedure for tracking the line uj{t) through the Grobner fan detecting 
when uj{t) leaves a cone. The salient point of the generic Grobner walk is that this calculation 
can be done formally by choosing certain generically perturbed uq and tq given by -<i and -<2. 

Here are the steps of the usual Grobner walk algorithm with the modified lifting step. 
Recall that a marked polynomial is a polynomial with a distinguished term, which is the 
initial term with respect to a term order -<. For a marked polynomial /, d{f) is defined in the 
natural way (cf. the definition of d^{f) in ^2.3|) . A marked Grobner basis over a term order 
-< is a Grobner basis over -< with all initial terms (with respect to -<) marked. For a marked 
Grobner basis we let d{G) = Uf(zcd{f). 

INPUT: Marked reduced Grobner basis G for / over a term order -<i, a term order -<2 along 
with ujq G G^-^{I) and Tq G G^^{I). 

OUTPUT: Reduced Grobner basis for I over -<2. 

(i) t = — oo. 

(ii) Compute_last_t. If t = oo output G and halt. 

(iii) Compute generators m^{G) = {in^((7) \ g ^ G} for in^(J) as 



veSg 

where Sg = {v e snpp{g) \ {u} \ t„„^ = t} and Qii^OC IS the marked term oi g ^ G. 
(iv) Compute reduced Grobner basis H for in^(/) over -<2 and mark H according to -<2. 



Use marking of H to mark H'. 

(vi) Autoreduce H' and put G = H' . 

(vii) Repeat from (fiij). 
Compute_last_t: 

1. Let V := {v E d{G) \ {uJo,v) > and (ro,f) < and t < where 

(cuo, v) 




Let 



H' = {f- f I / G H}. 



2. li V = 0, put t = oo and return. 

3. Let t := min{t„|f G V} and return. 
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Figure 1: C^i,^2 -<i=degrevlex and -<2=lex 

4 The generic Grobner walk 

In this section we show how certain generic choices of ujq and tq from ^lead to the generic 
Grobner walk algorithm. The crucial point is that step Q of the procedure Compute_Iast_t 
can be carried out formally using uq and tq from well defined perturbations given the term 
orders -<i and ^2- 

For an ideal I C R we let d{I) C Q" denote the union of d^{G), where G runs through 
the finitely many reduced Grobner bases for /. Let -<i and -<2 be two term orders given by 
Q-bases u = {ui, . . . ,a;„) and r = (ri, . . . ,r„) of Q" respectively. Observe that cu^ and 
are in the interior of the Grobner cones C^^(/) and C^^{I) respectively for sufficiently small 
positive rj. This follows from Lemma [2.1 1 Now define 

C^,,^^ = {w G M" I ^1 V and v -<2 0}. 

Here -<i, -<2 are extended to group orders on using uj and r. 

Example 4.1 Suppose that -<i is degree (reverse) lexicographic order and -<2 lexicographic 
order with y -<i^2 x. Then choosing u = ((1,1), (0, —1)) and r = ((1, 0), (0, 1)), we get -<i w 
imples (0,0) <;ex ("^i + ^'2, —^2) cind v -<2 implies (f 1,1^2) <iex (0,0). Intersecting the regions 
yielded gives (see Figure 1) 

C^,,^, = {(x, y) G I X + y > 0, X < 0}. 
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To fully understand the choice of 6 and e in the following we encourage the reader to 
compare with the computations in (*) and (**) below. Define 

Mt- = {{Ti, u)v I 2 = 1, . . . , n; M, t> e d{I)}. 

Corollary 12.21 shows that there exists sufficiently small positive 5 such that 

u^iv {ujs.u) < {ujs,v) (1) 

for u,v E Mt-. Suppose that 6 satisfies (H)). Now put 

Ns = {{ujs,u)v \u,v e d{I)}. 
Again by Corollary 12.21 we know that there exists sufficiently small positiye e such that 

u -<2V {t„ u) < (r„ v) (2) 

for u,v G A''^. Suppose now that we pick 6 according to ^ and e according to If 
V G d{I) n C^i,-<2 we put 



{us,v)-{t„v) ^ (r„ v) 

{ujs, v) 

If u,v E d{I) n C^i,-;2 then {us.u), {ujs,v) > and 

tu ^ tv 



{us, u) {lus, v) 

(Te, {us, V)u) < {t„ {us, u)v) 

{uJs,v)u ~<2 {ujs,u)v 
To evaluate -<2 above we see that 

{n, {uJs,v)u) < {Ti, {ljs,u)v) (**) 

{us, {Ti, u)v) < {us, {Ti, V)u) 
{Ti,u)v ^1 {Ti,v)u 

for i = 1, . . . ,n. Let T denote the matrix whose rows are ri, . . . , r^. By choosing 6 and e 
generically as above it follows that 

tu < ty <^=^ Tuv*' ^1 Tt>M* 

where -<i above refers to the lexicographic extension of on Z'^ to x ■ ■ ■ x Z". Here, Tuv* 
and TwM* are n x n matrices and we need to compare their rows. Notice that the comparison 
between tu and t„ does not involve 6 and e but only the term orders -<i and -<2- This leads 
us to define the facet preorder -< by 

u ~< V <^=^ tu < ty <^=^ Tuv* ^1 Tvu* (3) 

for u,v E d{I) n C^^^^2- 
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Example 4.2 Continuing the setup in Example \4 1\ if u = (mi,M2) and v = {vi,V2), then 




and the facet preorder is given by 

U ^ V <^=^ {ui'V ^1 Viu) V {{uiV = Viu) A {U2V -<i V2U)). 

If tu = ty then Tuv'^ = Tvu^ and * = vu^ since T is an invertible matrix. The identity 
uv^ = vu^ imphes that u and v are colhnear. Since u and v he in the same half space, m is a 
positive multiple of v. 

This has the nice consequence that the line uj{t) between ujs and intersects the cones 
in the Grobner fan in dimension > n — 1. Consider the typical situation, where v & C = 
C{v^vi^ . . . ,Vm) is chosen to minimize t^ as in the Grobner walk. Then uj(ty) is on a proper 
face F of C^. Since ty = t^ implies that m is a positive multiple of f for m G {f i, . . . , Vm}, we 
conclude that dimF = n — 1 i.e. F is a facet. 

The facet preorder -< defined in (jS)) may be inserted in the classical Grobner walk algorithm 
giving the generic Grobner walk algorithm completely removing the numerical dependence on 
the line u{t). Below, — oo(oo) denotes a vector strictly smaller (larger) than the vectors in 

INPUT: Marked reduced Grobner basis G for / over a term order -<i and a term order -<2 
(the facet preorder -< is given as in using and ^2)- 

OUTPUT: Reduced Grobner basis for I over -<2. 

(i) w = —00. 

(ii) Compute_last_w. If w = 00 output G and halt. 

(iii) Compute generators m^{G) = {m^{g) \ g G G} for in^(J) as 

V€Sg 

where Sg = {v & supp(5') \ {u} \ u — v^w,w^u — v} and a.^x^ is the marked term of 
geG. 

(iv) Compute reduced Grobner basis H for in^(/) over -<2 and mark H according to -<2- 

(v) Let 

H' = {f- f I / G H}. 

Use marking of H to mark H'. 

(vi) Autoreduce H' and put G = H'. 

(vii) Repeat from (jHl). 
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Compute_last_w: 

1. Let V ■.= {ve d{G) n C^,,^, | w -< v}. 

2. If = 0, put w = CO and return. 

3. Let w := min^{w|f G V} and return. 

4.1 Variations on the generic Grobner walk 

Several variations on the generic Grobner walk are possible. In many cases generators for 
an ideal are given which form a natural Grobner basis with respect to a specific weight vector. 
This happens for example in implicitization problems with polynomials yi — /i, • • • , — fm, 
where /, are polynomials in xi, . . . , x„ for i = 1, . . . , m. These polynomials form a Grobner 
basis over a vector uj assigning zero weights to xi, . . . , x„ and positive weights to . . . , Um- 
In this case one only needs to work with u and perturbations of the target vector. One 
may also truncate the facet preorder -< (to get a face preorder) using only parts (cui, . . . ,ujp) 
and (ti, . . . , Tg) of the Q-bases u and r. This leads to an analogue of the perturbation degree 
{p, g)-walk defined in j2] . 

5 An introductory example 

We illustrate the generic Grobner walk by a detailed example in the two dimensional case. For 
a given polynomial / e we let Lc^f) = f — f'^, where G is a marked Grobner basis (markings 
are underlined). Let -<i denote degree (reverse) lexicographic order and -<2 lexicographic order 
with y -<i^2 X. The facet preorder -< is given as in Example 14.21 Consider the ideal 

I = {x^ - y^, -y^ -x) C Q[x, y]. 

Initially we put 

G = {y^ - x^, - - x}, 

where the initial terms over -<i are marked. Clearly G is the reduced Grobner basis for I over 
-<i. The Grobner cone is given by 

C^,(/)=C({(-2,3)}U{(3,-2)})VnM|o- 

In this case (3, —2) ^ C'^i,-<2 and V = {(—2,3)}. So the first facet ideal is {y^ — x'^,x^). The 
reduced Grobner basis for this ideal over -<2 is {gf_ — y^, xy^, y^} and the lifting step is given 
by 

Lg{x' - 2/=^) =x'-y' 

Lcixy^) =xy^ -y"^ -X 
Lciy') = y' - xy^ - x\ 
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Our new marked reduced Grobner basis is 

G = {xi- ,x]^ - y"^ - x.i^ - xy"^ - y^}. 

Since w = (—2,3) ~< (—1,4) it follows that V = {(—1,4)} and the next facet ideal is 
(x^, xy^, — xy^) with reduced Grobner basis {x^, xy^ — y^, y'^} over -<2. Since 

Lg{x^) =x^ -y^ 
Laixy^ - 2/^) = xy^ - y"" + y"^ 

Lci.y'^) =y'^ -y"^ -y^ - X 

our new marked reduced Grobner basis is 

G = {^-y^xy^-y"^ + y\i-y^ -y^ -x}. 

Since w = (—1, 4) -< (—1, 7) we get V = {(—1,7)} and the next facet ideal is (x^, xy^, y'^ — x) 
with reduced Grobner basis {y^,x — y'^) over ^2- Here 

La{y') =y'- 2y^ -y' + y' 
Lg{x - y'^) = X - y'^ + y"^ + 

The new marked reduced Grobner basis is 

G = {y^-2y'-y' + y\x-y' + y^ + y'}. 

Since \^ = in this case, the generic Grobner walk halts and G is the reduced Grobner basis 
for / over -<2. 

6 Computational experience for lattice ideals 

In this section we report briefly on computations using the implementation GLATWALK ^3] 
of the Buchberger algorithm and generic Grobner walk for lattice ideals. Not surprisingly 
the walk performs best when initial and target vectors are close. An ideal situtation where 
this arises seems to come from a special case of feasibility in integer programming. Consider 
natural numbers Oi, . . . , a„ G N. Given 6 G N decide if the equation 

xitti H h x„a„ = b (4) 

has a solution xi, . . . ,Xn G N and find it if so. Adjoining the extra variable t we seek to 
minimize t subject to 

t + Xifli + ■ ■ ■ + XnCin = b (5) 

and t,xi, . . . ,Xn > 0. We denote this integer programming problem IPA^rib), where A is the 
1 X (n + 1) matrix (1 ai ... a„) and r = (1, 0, . . . , 0) G N""*"^. This problem has a trivial 
feasible solution: t = b,Xi = ■ ■ ■ = Xn = 0- Now we may apply standard algebraic techniques 
in integer programming (cf. jTj and ^HI) and form the toric ideal 

Ia = (xi -r\...,x„-r") C Q[t,xi,...,x„]. (6) 



14 



A Grobner basis G-r for I a with respect to r is a test set for the integer programming problems 
IPa,t{^)i where h varies and an optimal solution to (0) is the exponent of the normal form of 
with respect Gr thereby solving (jij). 
It is important to observe that the generating set for I a in (jHI) already is a Grobner basis 
G(j for I A with respect to the vector a = (—1, 0, . . . , 0). In the following section we report on 
computational results in computing Gr using the generic walk to go from cr to r compared with 
a direct computation with Buchberger's algorithm. We use the programs walk and gbasis of 
the program package GLATWALK. 

6.1 Comparison with Buchberger's algorithm 

To walk from cr to r we break ties with the reverse lexicographic order < given by t < Xi < 
■ ■ ■ < Xn i.e. we walk from the initial term order to the target term order <^. The names 
of the computational examples in the following table refer to specific numbers ai, . . . , a„ as in 
§ni They can be found in [Tj. The timings below are in seconds and the computations were 
carried out on a 1.6 GHz Pentium mobile with 1MB L2 Cache. 



EXAMPLE 


walk 


gbasis 


\Ga\ 


1 1 


cuwwl 


1.1 


17.7 


5 


7343 


cuww2 


11.4 


2.4 


6 


2472 


cuww3 


24.4 


9.5 


6 


4888 


cuww4 


1.2 


21.3 


7 


7937 


cuww5 


7.9 


1.3 


8 


1724 


probl 


0.1 


0.1 


8 


410 


prob2 


0.0 


0.0 


8 


142 


probS 


0.1 


0.1 


8 


425 


prob4 


0.1 


0.2 


8 


757 


prob5 


0.2 


0.1 


8 


516 


prob6 


0.1 


0.5 


10 


1035 


prob7 


0.1 


0.1 


10 


461 


probS 


0.2 


0.1 


10 


558 


prob9 


0.0 


0.0 


10 


270 


problO 


0.6 


2.5 


10 


2416 



In the problems cuwwl , cuww4 and problO the initial and target vectors are separated by 
less than 10 Grobner cones in the Grobner fan . This leads to surprisingly fast computation 
of relatively large Grobner bases. It would be interesting to further explore the efficiency of 
the generic Grobner walk in solving Frobenius problems more general than the Aardal-Lenstra 
knapsack problems. 
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7 Concluding Remarks 

The strength of the generic walk is that it is a completely deterministic algorithm avoiding 
the inherent instability of explicit numerical computation. 

The most recent version of the program 4ti2 (see jT2j) employs a projection algorithm for 
computing Grobner bases of toric ideals before using the usual Buchberger algorithm. This 
addition turns out to be a crucial optimization. It would be interesting to use a similar pro- 
jection algorithm before using the walk. Our preliminary experiments indicate that the walk 
is a little faster than the projection algorithm in the cases where it is significantly faster than 
the Buchberger algorithm. The current available implementations of Buchberger's algorithm 
like CoCoa, Singular, Macaulay2 etc. are slowed down by a significant factor compared to 
specialized integer vector implementations for lattice ideals. 

While the generic Grobner walk is presented here as a technique to compute a Grobner 
basis efficiently, one can use it for walking the entire Grobner fan systematically. In fact, a 
recent paper jTHj presents an algorithm based on both the generic Grobner walk and the reverse 
search technique to list all Grobner bases of a general polynomial ideal. In short, it reverses 
the generic Grobner walk in all possible ways from the lexicographic basis to reach all other 
bases. Obviously, such an exhaustive search requires an enormous amount of computational 
effort, and the symbolic perturbation turns out to be essential for this purpose. 

References 

[1] K. Aardal, A. Lenstra. Hard equality constrained integer knapsacks. Mathematics of 
Operations Research 29 (2004), 724-738. 

[2] B. Amrhein, O. Gloor, and W. Kiichhn. On the Walk. Theoret. Comput. Sci. 187 (1997), 
179-202. 

[3] B. Buchberger. An Algorithm for Finding a Basis for the Residue Class Ring of a Zero- 
Dimensional Polynomial Ideal (German). PhD-thesis (1965). 

[4] V. Chvatal. Linear Programming. W.H. Freeman and Company, 1983. 

[5] S. Collart, M. Kalkbrenner, and D. Mall. Converting bases with the Grobner walk. J. 
Symbohc Comp. 6 (1997), 209-217. 

[6] D. Cox, J. Little, D. O'Shea. Using Algebraic Geometry, 2. edition. Springer Verlag, New 
York 2005. 

[7] P. Conti, C. Traverso. Buchberger Algorithm and Integer Programming. Apphed Algebra, 
Algebraic Algorithms and Error-Correcting Codes, H. F. Mattson, T. Mora, T. R. N. 
Rao(eds), Lecture Notes in Computer Science 539, Springer- Verlag, (1991) 130-139. 

[8] M. de Berg, M. van Kreveld, M. Overmars, and O. Schwarzkopf. Computational Ge- 
ometry: Algorithms and Applications. Springer- Verlag, Berlin, Germany, 2nd edition, 
2000. 



16 



[9] J. Faugere, P. Gianni, D. Lazard and T. Mora. Efficient computation of zero- dimensional 
ideal Grobner bases by change of ordering. J. Symbolic Comp. 16 (1993), 329-344. 



[10 
[11 

[12: 

[13 
[14 

[15 
[16 
[17 
[18 
[19 
[20 



K. Fukuda, A.Jensen and R.Thomas. Computing Grobner fans, [math. AC/0509544 



T. Granlund. GNU Multiple Precision Arithmetic Library 4.1.4. September 2004. Avail- 
able at http : //swox . coin/gmp/| 

R. Hemmecke, R. Hemmecke, and P. Malkin. 4ti2 Version 1.2 — Computation of Hilbert 
bases. Graver bases, toric Grobner bases, and more. September 2005. Available at 
www.4ti2.de. 

N. Lauritzen, Truncated Grobner fans and lattice ideals, math.AG/0509247| 

N. Lauritzen, GLATWALK - software for experimentation with lattice ideals. August 2005. 
Available at [http : //home . imf . au . dk/niels/GL ATWALK/] 

T. Mora and L. Robbiano. The Grobner fan of an ideal. J. Symbolic Comp. 6 (1988), 
183-208. 

L. Robbiano. Term orderings on the polynomial ring. Proceedings of EUROCAL 85. 
Springer Lec. Notes Comp. Sci. 204 (1985), 513-517. 

B. Sturmfels. Grobner Bases and Convex Polytopes. University Lecture Series 8, Amer. 
Math. Soc, Providence, RI, 1996. 

R. Thomas. Algebraic methods in integer programming. Encyclopedia of Optimization 
(eds. C. Floudas and P. Pardalos). Kluwer Academic Publishers, Dordrecht, 2001. 

Q. Tran. A fast algorithm for Grobner basis conversion and its applications. J. Symbolic 
Comp. 30 (2000), 451-467. 

T. Yan. The geobucket data structure for polynomials. J. Symbolic Comp. 25 (1998), 
285-293 



